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H ' Abstract 



The neutron long counter NERO was built at the National Superconducting 
Cyclotron Laboratory (NSCL), Michigan State University, for measuring /3- 
delayed neutron-emission probabilities. The detector was designed to work 
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in conjunction with a /3-decay implantation station, so that (3 decays and 
/3-delayed neutrons emitted from implanted nuclei can be measured simulta- 
neously. The high efficiency of about 40%, for the range of energies of interest, 
along with the small background, are crucial for measuring /3-delayed neu- 
tron emission branchings for neutron-rich r-process nuclei produced as low 
intensity fragmentation beams in in-ffight separator facilities. 

Keywords: Large neutron counter, /3-delayed neutron emission, 
Astrophysical r-process. Neutron detection efficiency. Neutron background 
PACS: 28.20.-v, 28.20.Gd, 29.40.-n, 29.40. Cs, 23.40.-s, 25.40.Ny 



1. Introduction 

The emission of /3-delayed neutrons by neutron-rich nuclei significantly 
influences [l| the nucleosynthesis of heavy elements in the rapid (r-) neutron- 
capture process [2, ISl]. This decay mechanism competes with the /3 decay of 
r-process nuclei towards the valley of stability and serves as an additional 
source of neutrons in late stages of the r-process [^. Measurements of f3- 
delayed neutron emission probabilities {Pn) are needed for reliable r-process 
model calculations, and to test the astrophysical assumptions in various r- 
process models by comparing their final abundance predictions with obser- 
vations. 

From a nuclear-structure point of view, the P„ value provides model con- 
straints at low beam intensities where 7-spectroscopy is difficult. The P„ 
value probes /3-decay strength at excitation energies slightly above the neu- 
tron threshold. It therefore provides nuclear structure information comple- 
mentary to /3-decay, which often favors low energy /3-decay strength owing 
to the larger phase space (see for example p, |6[). 

The experimental determination of Pn requires the measurement of /5- 
delayed neutrons in coincidence with the /3 particles emitted from the nu- 
cleus of interest. This is particularly challenging for nuclei near or at the 
r-process path due to their very low production rates and the relatively short 
half-lives — of the order of 10-100 milliseconds. Experiments performed at 
ISOL-type facilities have successfully exploited the use of neutron long coun- 
ters (NLC) [7l| to measure Revalues of neutron-rich nuclei (see, for instance, 
the compilations of Refs. [Sl, loj). NLCs generally consist of a series of gas pro- 
portional counters embedded into a moderator block used to thermalize the 
neutrons prior to their detection. Performance requirements include a high 



detection efficiency for neutron energies ranging from a few keV to ~1 MeV. 
Because the detector does not measure the energy of individual neutrons, 
variations of the efficiency as a function of energy have to be minimized as 
much as possible as they otherwise can translate into uncertainties in the 
measured Pn- Our goal was to keep detector induced systematic uncertain- 
ties well below the 10% level. Measurements at that level of accuracy are 
a dramatic improvement over theoretical predictions, and ensure that other 
uncertainties dominate astrophysical and nuclear structure models. With 
systematic errors at that level, statistical errors will tend to dominate in 
practice, as the most interesting isotopes will typically be produced at rather 
low rates. 

We report here the development of NERO, a new NLC at National Super- 
conducting Cyclotron Laboratory (NSCL) suitable for use with fast radioac- 
tive beams produced by in-flight fragmentation. This technique provides 
exotic beams without some of the limitations induced by chemistry-based 
target-extraction techniques. The short time required to transport, sepa- 
rate, and identify the produced fragments, typically less than few hundred 
ns, makes it possible to study the very short-lived nuclei in the r-process. The 
fragments of interest are implanted in an active catcher that is part of the 
NSCL Beta Counting System (BCS). Implantation of a fragment and emit- 
ted /3 particles are detected event-by-event. The correlation of decays with 
a previously implanted nucleus requires large area highly pixelated catchers, 
typically double-sided silicon strip detectors (DSSDs). The challenge in the 
design of NERO was to include a large cylindrical cavity capable of accom- 
modating such a system, while still fulfllling the performance requirements 
for the detection efficiency. The final design was inspired by existing NLC 
detectors such as the Mainz Neutron Detector 10| . 



2. Technical aspects 

2.1. Design 

The detector system consists of a 60x60x80 cm^ polyethylene matrix 
(density 0.93(1) g/cm^) with its long symmetry axis aligned with the beam. 
Along the beam axis, the matrix has a cylindrical cavity with a diameter of 
22.8 cm to accommodate the BCS (see Fig. [H left). 

NERO includes three different types of cylindrical proportional counters 
manufactured by Reuter-Stokes: filled with ^He (models RS-P4-08 14-207 and 
RS-P4-0810-104), and filled with BF3 (model RS-Pl-1620-205) (see Tab. [D 



for details). Sixty of these detectors are arranged in three concentric rings 
around the central symmetry axis, allowing for a nearly Air solid angle cov- 
erage around the implantation detector (see Fig. [1], right). The optimum 
detector configuration was found using the MCNP code [ll| to calculate the 
neutron-detection efficiency for different geometries, moderating materials, 
and number and arrangement of various types of proportional counters. In- 
teractions of neutrons with the different detector materials were calculated, 
using the ENDF/B-VI 1^ cross-sections in the energy range 10~^ eV to 
20 MeV. The influence of different environments such as laboratory floor and 
wall configurations were investigated but were found to be negligible. Ac- 
cording to these calculations, most of the neutrons emitted from the center 
of NERO are detected in the innermost ring. Therefore, sixteen of the more 
compact and efficient ^He gas-filled proportional counters are mounted in 
the innermost ring at a radius of 13.6 cm. For the middle and outer rings 
at radii of 19.2 cm and 24.8 cm we use twenty and twenty-four BF3 propor- 
tional counters, respectively. The BF3 counters are longer allowing one to 
cover more solid angle, and their efficiency just compensates the decreasing 
efficiency of the inner ring with increasing neutron energy. 
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Figure 1: Schematic drawings of the NERO detector. Left: Side view showing the BCS 
chamber located inside of NERO with the DSSD at the central position. Right: backside 
showing the cylindrical cavity to house the BCS and the three concentric rings of gas-filled 
proportional counters. The labels A, B, C and D designate the four quadrants. 



To facilitate transportation and assembly, the polyethylene block is di- 
vided into an upper and lower half, and each half is subdivided in six equal 
parts along the longest symmetry axis. The twelve pieces are held together 
with eight stainless steel bolts. 



Table 1: Technical specifications of the NERO gas- filled proportional counters, (a) and 
(b) refer to the ^He detector models RS-P4-0810-104 and RS-P4-0814-207, respectively. 



Detector Active 


Radius 


Nominal 


Gas 


High 


Length 




Pressure 


Composition 


Voltage 


(cm) 


(cm) 


(atm) 




(+V) 



^He (a) 


25.0(2) 


1.3(2) 


10.2 


100% ^He 


1350 


3He (b) 


35.6(2) 


1.3(2) 


4.0 


100% ^He 


1100 


BF3 


50.8(1) 


2.5(2) 


1.2 


>96% lOB 


600 



2.2. Electronics 

The NERO readout channels are grouped in 4 quadrants with 15 chan- 
nels each (4 ^He counters and 11 BF3 counters). Figure |2] shows a schematic 
diagram of the NERO electronics for one quadrant. The proportional coun- 
ters detect the charged particles produced in the exothermic neutron-capture 
reaction ^He(n,p) or ^^B{n,a), respectively. Their signals feed 16-channel 
preamplifiers built at NSCL using Cremat CR-IOID miniature charge-sensitive 
preamp chips. The pre-amplified signals are sent into four 16-channel shaper 
and discriminator modules, designed at Washington University, St. Louis, 



and manufactured by Pico Systems [13|. These modules integrate indepen- 



dent shaping and discriminating circuits sharing the same input. Shaping 
times and pole- zero cancelation are adjusted for each channel by properly se- 
lecting the capacitances. The gain and threshold levels of the shaper/discriminator 
are adjusted via computer control over the CAMAC bus. 

The logic signals from the discriminator are recorded in scalers and in 
a 64-channel multi-hit (VME) TDC that is common for all quadrants. The 
TDC was programmed to work in start-gate mode, in which a gate signal, 
generated by a /3 decay detected in the BCS, enables the module to accept 
multiple stop signals in each channel from any of the sixty gas counters. The 
duration of this gate (r=200 /is) was chosen to account for the time needed 
to moderate and detect the neutrons (see Sec. 13. ip . The P„ value of a given 
nucleus is extracted from the number of stops-signals registered in the TDC 
(i.e., neutrons correlated with /3 decays) relative to the number of (3 decays 
detected in the BCS. 

The shaper outputs are connected to 32-channel (VME) ADC cards. The 



pulse height spectra recorded by the ADCs are used to set the thresholds 
and gains of the shaper/discriminator units, and to monitor any background 
or gain variation during the course of an experiment. Figure [3] shows typical 
ADC spectra for ^He and BF3 gas counters recorded under different condi- 
tions. The spectra show the typical wall-effect: The location of the peak 
at high amplitudes marks the Q value of the neutron-capture reaction, i.e., 
^B.e{n,p)t and ^°B(n,Q;)^Li for the ^He and BF3 gas counters, respectively; 
the plateau or low energy tail at low amplitudes arises from events where 
reaction products hit the detector wall preventing the complete deposition of 
their energy in the counter gas. Thresholds are set below these low amplitude 
events and just above the tail of the prominent low energy peak generated 
by electronic noise and background 7 radiation. Note that unlike the TDC, 
the ADCs only register one neutron per ECS trigger. 
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Figure 2: NERO electronic diagram (see text for details), 
quadrant is shown (GG stands for Gate Generator). 
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Figure 3: NERO ADC spectra for one of the "^He (left panels) and BF3 (right panel) gas 
counters. The top panels were recorded during 5 minutes, using a ^'^^Cf neutron source 
at the center of NERO. The middle panels correspond to a one-hour measurement with 
nuclei produced in the reaction ^^^Xe (150 MeV/u)-|-Be that were implanted in the BCS 
and included /3 delayed neutron emitters. The data displayed in the bottom panels were 
recorded during a 12-hour background measurement without beam on target using NERO 
as trigger. 



3. Detector performance 

3.1. Moderation time 

In order to determine the optimal duration of the TDC gate r, we mea- 
sured the time needed by the neutrons to slow down in the polyethylene 



moderator before their detection (r„). A ^^^Cf source was located at the 
center of NERO facing a Nal scintillator at a distance of 5 cm. Neutrons 
and 7 rays were emitted in coincidence from the fragments produced in the 
spontaneous fission of ^^^Cf. The scintillator was used to detect the 7 rays, 
which provided the external trigger of the NERO electronics (replacing the 
BCS trigger shown in Fig. [2]). The time difference between the detection of 
a 7 ray and the subsequent moderated neutrons recorded in the multi-hit 
TDC provided r„. The distribution of r^ is shown in Fig. H] (left) for the 
innermost, intermediate, and outer rings of proportional counters, as well as 
for the whole detector. 
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Figure 4: (Color online). Left: Measured background-subtracted moderation time dis- 
tributions for a ^^^Ci source for the innermost ring (squares), the intermediate ring (tri- 
angles), the outer ring (empty circles) and for the entire detector (solid circles). Right: 
Moderation-time distribution for each ring (empty circles), compared with results obtained 
with MCNP (solid lines). For easy comparison, the measured and calculated distributions 
for the first and second rings are scaled by 100 and 10, respectively. Note that the error 
bars of the experimental data are smaller than the symbol sizes. 



The largest differences between the three rings are found at the shortest 
times, when most of the neutrons emitted from the center of NERO reach 
the innermost ring. At late times, the neutrons are more uniformly dis- 
tributed over the whole moderator, and the three rings have similar detection 
rates. For each ring, the excellent agreement of the measured, background- 
subtracted moderation-time distributions with MCNP simulations is shown 
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in Fig. m Between 50 /is and 300 /is the time distributions can be approxi- 
mated with exponential functions. The corresponding measured and calcu- 
lated moderation time scales are 43 /is and 41 /is for the first ring, 51 /is and 
52 fis for the second ring, and 55 /is and 59 /is for the third ring, respectively. 
From Fig. |U we find that 94.3(1)% of the neutrons are detected within 
r„ <200 /is. The energies En of the neutrons emitted in the spontaneous 
fission of ^ssqj ^j^q typically described by a Maxwell-Boltzmann distribution 



function with an effective temperature of kT=lA2 MeV [1^, an average 
neutron energy of 2.1 MeV and a smooth tail at higher energies that extends 
up to 9 MeV. These are higher energies than typically expected for /3-delayed 
neutrons (see discussion at the end of Sec. I3.2.4p . Since the moderation time 
increases with neutron energy, we expect that more than 94.3(1)% of P- 
delayed neutrons are detected within 200 fis. We therefore chose a TDC gate 
of r=200 /is. 

3.2. Efficiency 

In order to characterize the NERO efficiency and its energy dependence, 
different types of measurements were performed using a ^^^Cf neutron source 
of known activity and neutrons produced in resonant and non-resonant reac- 
tions at the Institute for Structure and Nuclear Astrophysics (ISNAP) at the 
University of Notre Dame. The results obtained from these reactions were 
used to constrain the energy dependence of the NERO efficiency. 

3.2.1. Measurement of NERO efficiency with a '^^'^ Cf neutron source 

Before and after the calibration measurements performed at the Univer- 
sity of Notre Dame, the NERO efficiency was measured with a 1.251(5) /xCi 
^^^Cf calibration source with an active diameter of 5 mm (neutron branching 
11.6% and half-life 2.689 years). Additional contributions to the total neu- 
tron rate from contaminants were estimated. Besides ^^^Cf, there are small 
impurities of 249-25iQf ^.^^ 254^^ 249Qf g^j^j 25iQf ]iQyQ a negligibly small 
spontaneous fission branch, whereas the present amount of 254Qf ^^^ very 
small due to its short half-life. Consequently, only the ^soq^ ^^^^ "^^"^Ci iso- 
topes had to be considered. At the time of the measurement, 3.6% of the total 
activity of the source was due to ^^°Cf, whose contribution to the neutron 
activity was negligible due to its very low neutron branching of 0.296%. The 
alpha decay of 252Qf produces ^^'^Cm, which undergoes spontaneous fission 
accompanied with neutron emission, with a branching ratio of 8.39%. The 



very long half-life of this radioisotope (3.48x10^ years) made its contribution 
to the total neutron rate negligible. 

The number of detected neutrons was recorded with scalers and the multi- 
hit TDC described in Sec. 12.21 We verified that data processing dead time 
is negligible up to 50 KHz, well above the activity of the source, using a 
random pulser. Similarly, the 2 fis dead-time in the proportional counters 
was negligible. Taking the ratio of the number of neutrons recorded with 
NERO to the number of neutrons emitted by the source, calculated from 
the known source activity and neutron branching, we obtained a neutron 
detection efficiency of 31.7(2)%. This value also serves as a reference point 
to verify the NERO efficiency before, during, and after experiments. 

3.2.2. Measurement of NERO efficiency with resonant reactions 

The two resonant reactions used to st udy the NERO efficiency were 



^^C{a,nf^O [15|, [l6l ll7| and i^B(a,n)i^N [l8|, ll9|. The ISNAP KN Van 
de Graff accelerated a beam of a particles impinging onto ^'^C and ^^B tar- 
gets of 14(2) yug/cm^ and I2I2 /^g/cm^ thickness, respectively, located at the 
center of the NERO symmetry axis. The rate of incident a particles (/„) 
was monitored with an isolated electron-suppressed plate behind the target. 
A total of two resonances for the reaction ^^C(a,n)^^0, and one for the re- 
action ^^B(a,n)^^N were used (see Table [2]). Each resonance was completely 
mapped around its peak energy Er by detecting the number of neutrons 
as a function of a-beam energy E^. A linear background function was fit 
underneath the resonance curve and used to subtract non-resonant contri- 
butions and background neutrons. The efficiency of NERO was determined 
as the ratio of the number of detected neutrons to the number of neutrons 
Nn produced in the resonant reaction. Since the resonances considered here 
fulfill r^ << A£^ = Ei — Ef and F^ << Eji (see Table E]), one can use 
the thick-target narrow-resonance approximation and calculate Nn (see for 
example |20|) using 

^^ _ Ij7c^h\u-f)RNAP fdEy' 

f^AEn [dzj ■ ^'^ 

where [dE/dz) is the stopping power of the a particle in the target material, 



calculated with the SRIM-2000 code 21 in the center-of-mass frame; t is 



the duration of the measurement; A is the mole mass of the target; p is the 
target mass density; /i is the reduced mass of the system and (0^7) h is the 
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resonance strength. The resonance parameters used in the calculations were 



taken from Refs. [15|, ll6|, |17|, ll8|, |19| . They are summarized in Table [2J 

The results of the efficiency measurements are shown in Table H] for the 
three selected resonances, along with their corresponding evaluated average 
neutron energies in the laboratory frame {En) . The calculation of this latter 
quantity is described in Sec. 13.2.41 



Table 2: Properties of the three selected resonances: resonance energy in the laboratory 
frame -Eq, resonance width T n, strength (uj"f)fi, excitation energy E^^ and spin and parity 

Reaction Ea Tji (1^7) i? E^ -P 

(MeV) (keV) (eV) (keV) 



i3C(a,n)i60 


1.053 


1.5(2) 


11.9(6) 


7165 


5/2 


i3C(a,n)i60 


1.585 


<1 


10.8(5) 


7576(2) 


7/2 


iiB(a,n)i4N 


0.606 


2.5(5) X 10-3 


0.175(10) 


11436 


7/2 



Table 3: Validation of the thick-target, narrow-resonance approximation, dp is the target 
thickness and /S.E is the energy loss in the target at the resonance energy. 

Reaction E^ dp AE Tji 

(MeV) (mg/cm2) (keV) (keV) 



'^C{a,ny^O 


1.053 


0.014 


18 


1.5(2) 


i3C(a,n)i60 


1.585 


0.014 


16 


<1 


i^B(a,n)i^N 


0.606 


0.012 


18 


2.5(5)xl0-3 



3.2.3. Measurement of NERO efficiency with non-resonant reactions 

Additional measurements of the NERO efficiency were performed at IS- 
NAP using neutrons produced in the ^^V(p,n)^^Cr reaction |22| at three 
different energies. This reaction has been used in the past for neutron de- 



tector calibrations [16|, ll8|. Here, a proton beam was accelerated at the KN 
accelerator and impinged onto a ^^V target mounted in the center of NERO. 

11 



Three incident proton energies of 1.8 MeV, 2.14 MeV and 2.27 MeV were 
chosen from regions of the excitation function with no individual resonances, 
using three targets with a thickness of 32 /xg/cm^. 

To determine the number of ^^V(p,?T,)^^Cr reactions that have occured 
during a measurement one can take advantage of the fact that for every 
^^ V(p,ra)^^Cr reaction, a radioactive ^^Cr is created with a half-hfe of 27.7025(24) 
days. The electron-capture decay of ^^Cr is followed by the emission of sev- 
eral X-ra ys 23l| and a 320.0824(4) keV 7 ray J24| with a branching ratio of 



9.91(1)% [23[. The number of neutrons Nn produced in the reaction can then 
be simply determined from the activity of the target after irradiation. The 
number of 320.1 keV 7 rays emitted was measured offline in a lead-shielded 
setup, where the irradiated target was mounted in a plastic holder facing a 
HPGe 7 detector. Decay losses during the irradiation, transport, and offline 
counting were negligible. The HPGe efflciency at 320.1 keV was found to 
be 0.76(4)% using a ^^^Ba calibration source, and by interpolating the effl- 
ciencies measured for the two 7 rays emitted at 302 keV and 356 keV. The 
deduced NERO neutron efflciencies are listed in Table Hlfor the three different 
proton energies. The systematic error is dominated by the 5% uncertainty 
in the activity of the 7-ray calibration source. 

3.2.4- Results and discussion 

In order to evaluate the energy-dependence of the efflciency, the en- 
ergy spectrum of the emitted neutrons needs to be known for each reaction 
used. These neutrons are, for our purposes, essentially mono-energetic in 
the center-of-mass frame. The center-of-mass energies En can be calculated 
from the known reaction Q-values. However, the corresponding laboratory 
frame neutron energy En depends on the center-of-mass polar angle 6 of the 
emitted neutron with respect to the beam axis. This leads to a broadening 
of the neutron energy distribution, with an average neutron energy (En) of 

.+1 
(En) = I E^^{x)Wix)dx, (2) 

where x = cos 6. The angular-correlation functions W{x) for each of the three 
resonances used here were calculated as described in [25;] and are shown in 
Fig. O In the case of the ^^B(q;, n)^^N reaction, the coupling of the neutron 
spin with the ground-state of ^^N leads to two possible values of the final spin. 
Variations of the results obtained using the two possible W{x) functions were 
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Table 4: Laboratory frame projectile energy Eproj, average neutron energy (En), and 
width of the neutron energy distribution A_E„ with corresponding measured neutron de- 
tection efficiency e„. The last three columns list efficiencies for isotropic sources emitting 
neutrons at the energy (£"„). e„ "isotr" is derived from the experimental value e„ us- 
ing corrections calculated with MCNP, e„ "MCNP" is the calculated efficiency, and e„ 
"MCNP scaled " is the calculated efficiency scaled to obtain a best fit to the experimental 
data. 



Reaction E, 



pro] 



(En) 



AE„ 



(MeV) (MeV) (MeV) (%) 



^n Cn 6n 

isotr. MCNP MCNP 
scaled 

(%) (%) (%) 



"B(a,n) 


0.606 


0.56 


0.12 


33(2) 


38(2) 


42 


37 


^'^C{a,n) 


1.053 


2.8 


0.31 


24(1) 


30(1) 


29 


26 


^''C{a,n) 


1.585 


3.2 


0.41 


27(1) 


33(1) 


29 


24 


''V{p,n) 


1.80 


0.23 


0.014 


39(2) 


36(2) 


44 


39 


''V{p,n) 


2.14 


0.55 


0.024 


34(2) 


32(2) 


42 


37 


^^\{p,n) 


2.27 


0.68 


0.028 


34(2) 


34(2) 


41 


36 



included in the final uncertainty. For ^^V(p,n)^^Cr we have chosen energies 
where the excitation function shows non resonant behavior. This justifies 
the assumption of isotropic neutron emission in the center-of-mass frame, or, 
equivalent ly, W{x) = 1/2. 

The average neutron energies calculated in the laboratory frame (£"„) to- 
gether with the width of the energy distribution AEn are shown in Table SJfor 
each reaction, along with the corresponding measured efficiencies e„. In the 
case of ^^V(p,n)^^Cr AEn also includes a small contribution from the energy 
loss and straggling in the target. In order to determine the efficiencies for neu- 
trons emitted isotropically with a given energy from the DSSD catcher in the 
BCS, our measured e„ need to be corrected for the angular distribution and 



energy range of the neutrons (see e.g. Ref. f26| for the reaction ^^C{a,ny^O). 
This was done using MCNP simulations: First, the NERO efficiencies were 
calculated with MCNP at the energies (En) of Table lU assuming an isotropic 
mono-energetic neutron source located in the center of the detector. A sec- 
ond calculation was then performed, using the calculated laboratory frame 
angular and energy distributions of the neutrons. The MCNP-calculated 
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Figure 5: Left: Angular correlation W{x) calculated in the center-of-mass frame for the 
7.165 MeV, 5/2^ (solid line) and 7.576 MeV, 7/2~ (dashed line) resonances in the reaction 
^^C(a,ri)^^0. Right: Angular correlation W{x) calculated in the center-of-mass frame for 
Jj=l/2+ (sohd hue) and J/=3/2+ (dashed hue) of the 11.436 MeV, 7/2" resonance in 
the resonant reaction ^^B{a,n)^'^'N. 



anisotropic-to-isotropic efficiency ratios were then used as a correction factor 
to translate the measured efficiencies into efficiencies for isotropic emission 
at a single energy (En)- In Table HI we show the measured efficiencies for 
the different reactions (e„) and the corresponding corrected values (e„ isotr.). 
The strongest correction of about 15% arises mainly from the angular cor- 
relation W{x) in the resonant reactions. The isotropic efficiencies are shown 
in Fig. [6] for the whole detector (left), and for each ring separately (right). 

The experimentally-determined efficiencies covered a range of energies 
from about 0.2 MeV to 3 MeV. In order to extrapolate the results to lower 
energies, we combined the experimental values (e„ isotr.) with a MCNP 
calculation of the efficiency as a function of energy (see dotted line in Fig. |6l 
left). Despite a global absolute overestimation of about 5%, the calculated 
en{En) function follows very well the energy dependence obtained from the 
measured data. 

As shown in Fig. |6] (right), the efficiencies calculated independently for 
each ring follow reasonably well the energy trends of the experimental data. 
For all the reactions investigated, MCNP reproduces the efficiency of the 
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Figure 6: (Color online). Left: MCNP-calculated total efficiency as a function of the 
neutron energy En scaled to the experimental data (solid line) and unsealed (dotted line), 
compared with measurements for the reactions ^^B(Q,n) (solid circle), ^^C(a,n) (empty 
triangles), ^^V(p,n) (empty circles), and with the ^^^Cf neutron source (solid square). 
The energy width of the ^^^Cf measurement was calculated according to the shortest- 
interval criterion defined by W. Briichle for asymmetric distributions 27|. Right: MCNP- 
calculated efficiencies as a function of neutron energy _E„ for the innermost ring (solid line), 
intermediate ring (dashed line) and outer ring (dotted line), compared to the measured 
values for the first (solid circles), second (empty circles) and third (empty squares) rings. 
Note the different scales of the two figures. 



innermost ring, which is the most efficient of the three. The agreement is 
somewhat worse for the other rings at the lowest energies. In particular, the 
calculations overestimate the efficiency of the middle ring by about 3% at 
energies below 700 keV. We investigated the possibility that this discrepancy 
could be related to the type of detectors used. Several test measurements 
of the efficiency were performed with a ^^^Cf source, using one single ^He 
proportional counter placed in the first, second and third ring, using the ^^^Cf 
source. A comparison of the calculated efficiencies with the values measured 
under these conditions was consistent with the results shown in Fig. |6] (right). 
The overall 5% absolute overestimation of the efficiency by the calculations 
could not be attributed to the uncertainty in the polyethylene density. First 
of all, variations of the density modified the calculated efficiencies for the 
second and third ring in the opposite direction of the ffist ring. Secondly, 
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when the variations in density were hmited to the uncertainties provided by 
the suppher, no differences in the calculated results were observed. 

The good agreement of MCNP with experimental data observed in Fig. [6] 
(right) for the first NERO ring with a thinner moderator layer, and the small 
discrepancies found for the second and third rings with thicker moderator 
layers points to a limitation of MCNP to accurately calculate the scattering 
process of the neutrons in the moderator material. One possibility would be 
molecular vibrational and rotational excitation modes in the moderator ma- 
terial. Whereas this problem would be hardly observable in detectors with 
thin moderators (e.g. Ref. [28|), it would become more severe for thicker 
moderators. Interestingly, similar conclusions were drawn when comparing 
MCNP calculations with neutron-flux measurements performed with thick 



neutron detectors |29l . |30|- In order to compensate for these model defi- 
ciencies we scaled the calculated efficiencies for each ring independently to 
better match the experimental data (see scaled efficiency in Table H]). The 
new scaled efficiency-curve (solid line in Fig. [6l left) can thus be used to 
extrapolate the efficiency to energies below 200 keV. 

It is worth noting that the relevant neutron energy range for /3-delayed 
neutron emission in r-process nuclei is a few hundred keV. As an example, 
for the r-process nuclei around A ~ 100 — 130, spectroscopic studies of (3- 



delayed neutron emitters |31l . |32| . |33| showed that En is typically much lower 
than Qi3 — Sn- Neutron energies in the laboratory frame En were found to be 
199 keV for ^^Br, 450 keV for ^^Rb, and 579 keV for ^^'^l. This result was fur- 
ther supported by the measured average neutron energies of fission fragments 
from 235U {{En) =575 keV) and ^s^Pu {{En) =525 keV), where, in addition. 



very few neutrons were found at En^SOO keV |3ll . |34| . According to these 



authors, the reason for the "compressed" En spectra is the preferred popu- 



lation of the lowest excited states in the final nuclei [33| . Our experimental 
efficiency calibration therefore covers the most critical energy range, and the 
condition of an energy independent efficiency for /3-delayed neutrons is well 
fulfilled. As an example, e^ shown in Fig. |6] (left) shows a relative variation 
of about ±5% for energies below 800 keV. This variation will contribute to 
the final uncertainty of the measured P„. This uncertainty can be reduced if 
the neutron energies En can be constrained from experiment or theory. 

3. 3. Background 

One limitation for the measurement of -P„, particularly for very exotic 
nuclei, is the neutron background rate {Bn)- Its estimation requires to dis- 
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tinguish two different origins. First, there is the "intrinsic" background as- 
sociated with the electronics of the NERO detector and its sensitivity to the 
neutrons present in the environment (mainly cosmic rays). Secondly, dur- 
ing the course of an experiment, there are beam-induced neutrons produced 
by nuclear reactions unrelated to the /3-delayed neutron emission of inter- 
est. During experiments, neutron background rates measured with NERO 
in self-trigger mode can vary within about 5-10 s~^ depending on whether 
or not the beam is on target. As will be discussed later, the impact of these 
background rates is dramatically reduced when the neutrons are measured 
in coincidence with /3 decays. 

Analysis of the ring-counting ratios for background runs (self-trigger mode) 
and production runs (external trigger mode) support the idea of an external 
and a beam-induced background source. As shown in Fig. [7] (left), measure- 
ments performed with NERO in external trigger mode (i.e. from /3 decays in 
the BCS) with /^-delayed neutron emitters showed that the neutron count- 
ing rates were higher for the innermost ring and systematically decreased for 
the outer rings, in agreement with the results shown in Fig. [H] (right). On 
the other hand, background runs with beam off showed the opposite trend, 
with high rates in the outer ring, gradually decreasing for the inner ones 
(Fig. [TJ center). This result suggests that these runs were mainly affected by 
an external source of background neutrons, most probably related to cosmic 
rays. Finally, background runs with beam on target showed an intermediate 
situation that could be explained as arising from a combination of external 
and internal sources (Fig. [TJ right). Energy spectra obtained for background 
runs with the ADCs show the wall-effect shape expected for neutrons (see 
lower spectra in Fig. [3]). Electronic and 7-ray contributions are largely below 
the discriminator thresholds. 



4. Measurement of P, 



n 



The NERO detector, together with the BCS, has been eniployed in numer- 
ous r-process motivated experiments performed at NSCL p, la, |35|, l36|, |37 



The exotic nuclei of interest are implanted in a 40 x 40-pixel DSSD in the 
BCS. /3-decays are also detected in the DSSD and can be position-correlated 
to previously implanted ions during a maximum correlation rime tc- Pn values 
were determined using the number of neutrons Nj^n detected in coincidence 
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Figure 7: Ratio of neutrons detected with different NERO rings for three different runs: 
production with /3-delayed neutron emitters (left), background with beam off (center), 
and background with beam on (right). Histogram bin numbers 1, 2 and 3 correspond to 
ring ratios R2/R1, R3/R2, and R3/R1, where R1-3 are the innermost, intermediate and 
external rings. Statistical errors are negligible. 



with an implantation-correlated /3-decay event, according to the equation: 



Pn 



N, 



Br- 



N, 



enNp 



(3) 



where Bj^n is the number of background coincidences between /3-like events 
(including real /3 decays and background from the BCS) and neutrons, and iVg 
is the number of /3-decaying mother nuclei. Nf^pn is the number of detected 
/3-delayed neutrons from the daughter nuclei and needs to be subtracted from 
Npn- For the nuclear species analyzed in Refs. p, ly, |35|, |36|, |37| , /3-neutron co- 
incidences associated with descendant nuclei other than the ^-decay daughter 
were negligible. In this case, using the Batemann equations |38| , it is possible 
to write explicitly the value of Np^n as: 



where C is a constant given by: 



C 



A2 — Ai 



:i - Pn)C, 



1 - e-^^*^ - ^ (1 - e-^^*^ 
A2 



(4) 



(5) 



In this equation, P„„ is the neutron-emission probability of the daughter 
nucleus, and Ai and A2 are the decay constants of the mother and daughter 



nuclei, respectively. Inserting Eq. H] and Eq. E] into Eq. |31 and rearranging 
terms: 

_ Niin - -B/3n - C , , 

The value of N^ for a given nucleus is calculated as the product of the total 
number of implantations in the DSSD, and the /3-detection efficiency. 

The NERO background rate given in Sec. 13.31 is the "free" neutron back- 
ground rate without any coincidence requirements. In practice, however, P„ 
values are determined from neutrons measured in coincidence with /3 decays. 
The number of background /3-neutron coincidences for a given nucleus can be 
written as i?^„ = Bn{[i) + Bn{Bp); where Bn{f3) is the number of "free" back- 
ground neutrons in random coincidence with parent /3 decays, and Bn{Bp) 
is the number of "free" background neutrons in random coincidence with 
background /9-like events in the BCS. 

Bn{B/j) as a function of time and detector pixel can be reliably estimated 
from the /3-neutron coincidence rates outside of the correlation window of 
any ion implantations. Bn{Bp) for a specific parent nuclide can then be 
calculated by summing the specific backgrounds at the time and location of 
each individual ion implantation event. The high granularity of the DSSD 
detector greatly reduces this background. 

Bn{/3) can be calculated as the product of the number of parent (3 decays 
detected A^^ and the probability for at least one "free" background neutron 
to be detected in random coincidence with each parent /3 decay, i.e.: 

BniP) = iV, f; ^^e-^^^ = N,{1- e-^n, (7) 

fc=i 

where r is the TDC time window defined in Sec. 13. ![ and Rb is the "free" 
neutron background rate. For a typical value Rb —10 s~^, about 0.2% of 
the detected parent /3 decays are in random coincidence with a background 
neutron, setting the order of magnitude of the lowest P„ values that can be 
measured with NERO under these conditions. 

The Pn values and their errors obtained in various NSCL experiments can 
be found in Refs. p, |6|. The P„ values measured in these experiments agree 
well with perviously measured well established data. 
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5. Summary and conclusions 

The neutron detector NERO has been built at NSCL enabling the mea- 
surement of /3-delayed neutron emission probabilities of r-process nuclei with 
fast rare isotope beams. The specific design was motivated by the require- 
ment of achieving a high, energy-independent neutron detection efficiency 
up to ~ 1 MeV, and accommodating a large pixelated /3-counting system, 
necessary to perform measurements with fragmentation beams. MCNP sim- 
ulations were carried out during the design phase to find the optimum con- 
figuration. 

Studies of the detector efficiency at various neutron energies were per- 
formed with a 252Qf source, and with neutrons from a number of resonant 
and non-resonant reactions at ISNAP. MCNP calculations reproduce rea- 
sonably well the energy dependence of the detector efficiency. On the other 
hand, the MCNP calculations slightly overestimate the absolute efficiency of 
the second and third rings, when the neutrons traverse a larger volume of 
polyethylene. An overall scaling of the calculated efficiency to the measured 
data can be used to extrapolate the detector efficiency to smaller and larger 
neutron energies. The small energy dependence (of about 5%), for neutron 
energies below 800 keV, represents the main contribution from the efficiency 
correction to the total uncertainty of Pn- 

NERO is currently used with the BCS at NSCL, but can be used with 
other /3-decay stations at other rare isotope facilities. It will also be used to 
fully exploit the much higher production rates expected in new generation 
facilities like FRIB at NSCL, FAIR at GSI and RIBF at RIKEN. 
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